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ABSTRACT 

The late-type spiral galaxy NGC 6946 is a prime example of molecular gas dynamics 
driven by ‘bars within bars’. Here we use data from the BIMA SONG and HERA¬ 
CLES surveys to analyse the structure and stability of its molecular disc. Our radial 
profiles exhibit a clear transition at distance R ~ 1 kpc from the galaxy centre. In 
particular, the surface density profile breaks at R « 0.8 kpc and is well fitted by a 
double exponential distribution with scale lengths Ri « 200 pc and i? 2 ~ 3 kpc, while 
the ID velocity dispersion a decreases steeply in the central kpc and is approximately 
constant at larger radii. The fact that we derive and use the full radial profile of a 
rather than a constant value is perhaps the most novel feature of our stability anal¬ 
ysis. We show that the profile of the Q stability parameter traced by CO emission is 
remarkably flat and well above unity, while the characteristic instability wavelength 
exhibits clear signatures of the nuclear starburst and inner bar within bar. We also 
show that CO-dark molecular gas, stars and other factors can play a significant role 
in the stability scenario of NGC 6946. Our results provide strong evidence that grav¬ 
itational instability, radial inflow and disc heating have driven the formation of the 
inner structures and the dynamics of molecular gas in the central kpc. 

Key words: instabilities - ISM: kinematics and dynamics - galaxies: individual: 
NGC 6946 - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: structure. 


1 INTRODUCTION 

Non-axisymmetric galactic structures such as bars and spi¬ 
ral arms contain a significant fraction of the disc material, 
and thus have a prominent role in the evolution of their 
host galaxies (e.g., Lindblad 1960). Such structures are ef¬ 
ficient drivers for the redistribution of angular momentum, 
and induce strong flows and streaming motions in their hosts 
(e.g., Schwarz 1984). The interstellar gas is highly inhomo¬ 
geneous and prone to gravitational instability, thermal in¬ 
stability and turbulent compression (see, e.g., Mac Low & 
Klessen 2004). When large-scale flows converge, molecular 
clouds can form rapidly in regions dense enough to form H 2 
(e.g., Glover & Mac Low 2007; Tasker 2011). Galactic-scale 
flows also pump turbulent energy into the rotating gas disc, 
which in turn initiates local instabilities (e.g., Renaud et al. 
2013) and produces further heating of the disc (e.g., Martig 
et al. 2014). A thorough understanding of such multi-scale 
processes requires a detailed analysis of the structures that 
originate from them. Galaxies hosting a hierarchy of inter- 
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linked non-axisymmetric structures are thus especially suit¬ 
able for analysing disc instability and its effects on structure 
formation and evolution. 

The late-type spiral galaxy NGC 6946 is a prime exam¬ 
ple of such dynamics: it hosts an outer bar that drives the 
evolution of structure in the disc, from its prominent spiral 
arms to an inner bar within bar feeding a nuclear starburst 
(e.g., Elmegreen et al. 1992; Schinnerer et al. 2006; Fathi et 
al. 2007; Tsai et al. 2013). The gravitational instability of 
NGC 6946 has been investigated in several works, most often 
in the context of large galaxy surveys (e.g., Kennicutt 1989; 
Ferguson et al. 1998; Fuchs & von Linden 1998; Martin & 
Kennicutt 2001; Leroy et al. 2008; Romeo & Wiegert 2011; 
Hoffmann & Romeo 2012; Romeo & Falstad 2013). Fergu¬ 
son et al. (1998) showed how sensitive the results are to the 
gas ID velocity dispersion cr g : assuming that <r g = 6 kms -1 , 
the value proposed by Kennicutt (1989), NGC 6946 turns 
out to be unstable up to the edge of the optical disc, while 
using a radial profile of <r g derived from observations yields 
stability across the entire disc! Martin & Kennicutt (2001) 
pointed out that radial variation in cr g remains controver¬ 
sial because such measurements demand both high angular 
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resolution and high brightness sensitivity, requirements not 
met by most observations. Fortunately, recent CO and Hi 
galaxy surveys (BIMA SONG, HERACLES and THINGS) 
have provided high-quality measurements of gas kinemat¬ 
ics, which allow deriving reliable radial profiles of <r g (e.g., 
Caldu-Primo et al. 2013). 

Another source of concern in this stability context is the 
multi-component (gas+stars) nature of the disc. Unfortu¬ 
nately, there are no stellar spectra available for NGC 6946, 
i.e. we have no information about the stellar radial veloc¬ 
ity dispersion <r*. Romeo & Falstad (2013) used a model- 
based radial profile of <r*, originally proposed by Leroy 
et al. (2008), and analysed the stability of NGC 6946 to¬ 
gether with other spirals from The HI Nearby Galaxy Sur¬ 
vey (THINGS). They showed that NGC 6946 is unstable 
only within the central kpc, and that molecular gas (H 2 ) 
is the main driver of instability if oh 2 is within the range 
of values used in previous investigations, i.e. 6 kms -1 (e.g., 
Kennicutt 1989; Wilson et al. 2011) to 11 kms -1 (e.g, Leroy 
et al. 2008). They also found that if <th 2 ~ 6 kms -1 , then 
molecular gas controls the local stability level of the disc up 
to about 5 kpc from the centre. Stars dominate the value of 
the Q stability parameter at larger distances, while atomic 
gas has a negligible effect up to the optical radius. 

The multi-component stability analysis of Romeo & Fal¬ 
stad (2013) has one main limitation: it uses observationally 
motivated values of <th 2 , rather than a radial profile <th 2 (R) 
derived from observations. This is a serious limitation be¬ 
cause, as discussed above, molecular gas plays a primary 
role even beyond the central kpc, and up to a radius f? max 
that depends on 0 H 2 : -Rmax ~ 5kpc for <th 2 ~ 6 kms -1 ; 
the larger <th 2 , the smaller f? max . This illustrates how im¬ 
portant it would be to have a reliable radial profile of <th 2 
up to R ~ 5 kpc, and suggests that a one-component sta¬ 
bility analysis based on such a profile could be a proper 
starting point for discussing the roles that stars, the various 
gas phases, the various bar structures and other components 
play in the stability scenario of NGC 6946. 

In this paper, we analyse NGC 6946 focusing on the 
structure and stability of its molecular disc. For this pur¬ 
pose, we derive high-quality radial profiles of the surface 
density, ID velocity dispersion and epicyclic frequency up 
to R = 5 kpc, using data from two recent surveys: the 
BIMA Survey of Nearby Galaxies (BIMA SONG, Heifer 
et al. 2003), and the HERA CO-Line Extragalactic Survey 
(HERACLES, Leroy et al. 2009). We also analyse the link 
between structure formation and gravitational instability in 
the central kpc, and discuss the roles that CO-dark molecu¬ 
lar gas, stars and other factors play in the stability scenario 
of NGC 6946. 


2 NGC 6946 

NGC 6946 is a low-inclination grand-design barred spiral 
galaxy with an observed axis ratio of 0.87 (Carignan et al. 
1990; Zimmer et al. 2004; Boomsma 2007). We adopt a dis¬ 
tance of 5.9 Mpc (Leroy et al. 2009), so that l”corresponds 
to 28.6 pc. Numerous morphological and kinematic stud¬ 
ies by e.g. Elmegreen et al. (1992), Regan & Vogel (1995), 
Elmegreen et al. (1998), Kennicutt et al. (2003), Schinnerer 
et al. (2006) and Fathi et al. (2007) have revealed clear sig- 
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Figure 1. A composite BVI image of NGC 6946, with the North- 
East direction displayed in the top-right corner. The three ellipses 
show the outermost radius of the analysis presented here (5 kpc), 
and the approximate radii of the inner bar within bar (1 kpc 
and 200 pc). The ellipses take into account the inclination and 
position angle of the galaxy. Note how the dust lane near the 
centre in the North direction changes angle at the 1 kpc radius, 
i.e. approximately where the oILR of the outer bar is located. 
Image courtesy of Subaru Telescope (NAOJ) and Robert Gendler. 


natures of three main gravitational distortions: a large but 
weak bar-like structure with a projected ellipticity of 0.25 
and radius R « 4.5 / , a secondary bar with radius R « 60”, 
and a nuclear bar with ellipticity 0.4 and radius R « 8”. 
Fathi et al. (2007, 2009) found that the pattern speed of 
the weak outermost bar is lf p = 25 ± 6 kms -1 kpc -1 , 
and used this value to place the corotation radius between 
250” and 300”, the outer inner Lindblad resonance (oILR) 
radius between 30” and 60”, and the inner inner Lindblad 
resonance (ilLR) radius below 15”. In addition, direct mea¬ 
surements of the pattern speed of the secondary bar showed 
that this lies in the range 39 ± 8 kms -1 kpc -1 (Zimmer 
et al. 2004), 47^2 kms -1 kpc -1 (Fathi et al. 2007), and 
51 ± 3 kms -1 kpc -1 (Font et al. 2014). 

These dynamical and structural properties agree well 
with those found by Schinnerer et al. (2006, 2007), and con¬ 
verge towards a scenario where the outer bar drives the evo¬ 
lution of structure in NGC 6946 (see Fig. 1): the secondary 
bar is within the oILR of the outer bar, while the nuclear bar 
is within the ilLR of the outer bar and has approximately 
the same radial extent as the nuclear starburst (R « 5.5”; 
Schinnerer et al. 2006). Such structural interplay is indeed 
interesting from a dynamical point of view, and makes NGC 
6946 an ideal laboratory for a detailed analysis of its disc 
structure and stability. 
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3 DATA AND METHOD 
3.1 Molecular gas data 

In this work we analyse the molecular gas content and kine¬ 
matics of NGC 6946 as traced by the CO J( 1 —> 0) and 
CO J(2 —> 1) lines, and based on archival single-dish and 
interferometric observations from the BIMA SONG (Heifer 
et al. 2003) and HERACLES (Leroy et al. 2009) surveys, 
respectively. We have retrieved the reduced data cubes from 
the archives made available by these teams. All the details 
concerning data acquisition, data reduction and data qual¬ 
ity are presented in the two papers above (see also Leroy et 
al. 2008 for the HERACLES data). Most important for our 
analysis is that the spatial resolution of the CO J(1 —» 0) 
data is approximately 5" x 6”, sampled at l"/pix and with 
a channel width of 10 kms -1 . The CO J(2 —» 1) data have a 
spatial resolution of approximately 11", sampled at 2"/pix 
and with a channel width of 5.2 kms -1 . 

To derive the molecular gas surface density (E), line- 
of-sight velocity (Mos) and velocity dispersion (ct), we use 
two methods in parallel, namely calculating the moment 
maps and applying Gaussian fits to the individual spectra. 
We find a good agreement between the moment maps and 
the results of Gaussian fitting. We also find that our ampli¬ 
tude maps match the zeroth-moment maps presented by the 
BIMA SONG and HERACLES teams. One advantage of ap¬ 
plying Gaussian fits is that we can get better covering maps 
if we use a Hanning smoothing algorithm in the spectral di¬ 
rection before fitting the individual Gaussian profiles (e.g., 
Hernandez et al. 2005; Daigle et al. 2006). This procedure 
does not affect the line amplitude or shift. However, it arti¬ 
ficially broadens the individual spectra by Act w 4.4 kms -1 , 
which we then subtract quadratically from the derived veloc¬ 
ity dispersion maps. Finally, we apply a cleaning procedure 
by removing all the spectra for which the emission line am¬ 
plitude is smaller than twice the rms noise, ct is smaller than 
half the velocity channel, and the formal error in the derived 
ct is greater than 10 kms -1 . 

A careful inspection of the individual spectra reveals 
the presence of multiple components at different locations 
(mostly in the central few 100 pc). In view of the good agree¬ 
ment between the moment maps and the single-Gaussian 
fits, and in view of the complications involved in multiple- 
profile fitting schemes (Blasco-Herrera et al. 2010), we do not 
use profile decomposition methods. All subsequent analysis 
is thus based on single-Gaussian profile fitting. 


3.2 Radial profiles by means of robust statistics 

We derive the radial profiles of E and ct using robust statis¬ 
tics, which are especially useful when the data contain a 
significant fraction of outliers (see, e.g., Rousseeuw 1991; 
Muller 2000; Huber & Ronchetti 2009; Feigelson & Babu 
2012). This is indeed the case for the E and ct maps, where 
a non-negligible number of pixel values (e.g. those associated 
with Hu regions) deviate strongly from a normal distribu¬ 
tion. To derive the radial profiles of E and ct, we divide their 
maps into tilted rings, which are circular in the plane of the 
galaxy, and compute the median values of E and ct in each 
ring. We then estimate the uncertainty in these median val¬ 
ues via the median absolute deviation (MAD): 
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Figure 2. The rotation curve (top) and frequency diagram 
(bottom) of NGC 6946. Adopting a pattern speed of 25 ± 
6 kms -1 kpc -1 for the weak outer bar (Fathi et al. 2009), we 
recover resonance radii consistent witli the literature. The ro¬ 
tation curve is consistent with measurements based on different 
molecules, neutral atomic and ionized gas. The frequency diagram 
is based on smoothed curves so as to mimic the resolution of the 
curves presented in the literature (see Sect. 2). 

AAmed = 1.858 x MAD/n/JV , (1) 

MAD = median{|X; — X med |} , (2) 

where Xi are the individual measurements of E or ct, X med 
is their median value, and N is the number of resolution 
elements in each ring (i.e. the number of pixels in the ring 
divided by the number of pixels per resolution element). Eqs 
(1) and (2) are the robust counterparts of the formula tra¬ 
ditionally used for estimating the uncertainty in the mean: 
AA mean = SD / V~N, where SD denotes the standard devia¬ 
tion (see again Muller 2000). Indeed, the median and the 
median absolute deviation provide robust statistical esti¬ 
mates of the ‘central value’ and the ‘width’ of a data set, 
respectively, even when almost 50% of the data are outliers! 

Our numerous tests with varying ring radii and widths 
confirm that the high quality of the BIMA SONG and HER¬ 
ACLES data allow a derivation of E and ct at radii smaller 
than the synthesized beam size of the interferometric ob¬ 
servations. This is mainly thanks to the good sampling of 
the resolution element. The smallest reliable step is found 
to be 2", which corresponds to 57 pc. This is the step size 
adopted throughout our analysis. The inner 2.5 kpc of the 
E and ct profiles presented here are derived from the BIMA 
SONG data, while the profiles at larger radii are derived 
from the HERACLES data. The H 2 surface density is con¬ 
verted to physical units by adopting an Xco conversion fac¬ 
tor of 2 x 10 20 cm -2 (Kkms -1 ) -1 , consistent with Leroy et 
al. (2008) and Donovan Meyer et al. (2012). We also correct 
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for the contribution of Helium multiplying by a factor of 
1.36. To further match the BIMA SONG and HERACLES 
data, we use CO J(2 —> l)/CO J(1 —> 0) = 0.8, as in Leroy 
et al. (2008). 

In addition to the E and <r profiles, a stability analysis 
also requires the epicyclic frequency k at each radius, which 
we derive from the observed rotation curve. To calculate the 
rotation curve, we assume that circular rotation is the dom¬ 
inant kinematic feature, and that our measurements refer to 
positions on a single inclined disc. We then use the tilted ring 
method combined with the harmonic decomposition formal¬ 
ism (e.g., Schoenmakers et al. 1997; Wong et al. 2004; Fathi 
et al. 2005). Given that this procedure involves fitting sev¬ 
eral parameters at each radius of the observed velocity field 
(contrary to simply finding the median values of E and a), 
we cannot use the same initial step size. Hence we apply a 
larger radial step, and interpolate linearly to obtain the ro¬ 
tation curve at all radii where we have calculated the robust 
E and a values. Here the cut in the data is made at a radius 
of 0.7 kpc. Once the rotation curve Kot is calculated at each 
radius R, we derive the angular frequency f1 = V TO t/R and 
the epicyclic frequency k = \JRdQ 2 /dR + 4 Q 2 . Setting up 
the standard frequency diagram (see Fig. 2), we confirm the 
location of the inner Lindblad resonances of the weak outer 
bar discussed in Sect. 2. 


4 RESULTS 

4.1 Structure of the disc 

Fig. 3 illustrates that the radial profiles of the epicyclic fre¬ 
quency, k, ID velocity dispersion, cr, and surface density, E, 
break at R ~ 1 kpc. This tells us that the inner region of 
NGC 6946 is physically distinct from the rest of the disc, 
and that radial inflow (increase of S and k towards the cen¬ 
tre) and disc heating (increase of a towards the centre) have 
both played a primary role in shaping the central kpc. This 
is consistent with the fact that the dynamics of NGC 6946 
is driven by three bars, and that the strong secondary bar 
extends up to R ~ 1 kpc (see Sect. 2). Since bars are non- 
axisymmetric gravitational instabilities, they exert torques 
on the disc, which transport angular momentum and energy 
outwards, leading to the twoprocesses mentioned above: ra¬ 
dial inflow and disc heatingj (e.g., Zhang 1998; Griv et al. 
2002; Romeo et al. 2003, 2004; Fathi et al. 2008; Agertz et al. 
2009; Krumholz & Burkert 2010; Forbes et al. 2012, 2014). 
Besides, bars within bars make such processes highly effi¬ 
cient (e.g., Shlosman et al. 1989; Friedli & Martinet 1993; 
Rautiainen et al. 2002; Shlosman 2002). A dynamically cool 
inner disc and especially nuclear star formation make such 
structures long-lived, even in the absence of mode coupling 
(Du et al. 2015; Wozniak 2015). 


1 Disc heating is a natural consequence of radial inflow and is 
mediated by local gravitational instabilities. Although there are 
still open questions, the basic idea behind this process is simple, 
and is beautifully illustrated in sect. 7.1 of Kormendy & Kenni- 
cutt (2004). Radial inflow increases both E and re, but E ‘wins’ 
and the Toomre (1964) parameter Q = kct/itGT, decreases. As 
Q drops below a critical value of order unity, local gravitational 
instabilities set in and increase a , thus heating the disc. 
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Figure 3. Radial profiles of the epicyclic frequency (top), ID ve¬ 
locity dispersion (middle) and surface density (bottom) of molec¬ 
ular gas in NGC 6946. The thick lines are fits to these profiles over 
the radial range 1-5 kpc, extrapolated into the central kpc. The 
thin lines represent the approximate radial extent of the nuclear 
starburst and inner bar within bar. The bottom panel also shows 
a robust, median-based, exponential fit to S (R) for R < 0.8 kpc. 


Look now at the thick lines in Fig. 3. These are simple 
but accurate fits to the radial profiles of re, a and E for 
1 kpc < R < 5 kpc. 

(i) ttflt {R) is the fit originally proposed by Leroy et al. 
(2008); see their eqs (13) and (Bl), and table 4. 

(ii) crfit(-R) — 10 kms' 1 is the median of the ID velocity 
dispersion data. This value is comparable to those found 
by Walsh et al. (2002) and Caldu-Primo et al. (2013) for 
R> 2kpc. 

(iii) Efit(iT) = T, 2 e~ R / R2 , with E 2 ~ 66 M 0 pc~ 2 and 
i ?2 — 2.9 kpc, is a median-based fit to the surface density 
data. The disc scale length R 2 is 50% larger than that found 
by Leroy et al. (2008, 2009). Note, however, that they fit¬ 
ted an exponential function to the whole molecular disc, i.e. 
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they included the central kpc, which we have shown to be 
physically distinct from the rest of the disc. 

Extrapolating these fits back to R < 1 kpc allows us to pre¬ 
dict how the inner region of NGC 6946 differs dynamically 
from the rest of the disc. In particular, we have compelling 
evidence that the nuclear disc is also exponential. A robust, 
median-based, fit to the radial profile of E for R < 0.8 kpc 
yields: Efit(R) = Eie -fl / fll , with Ei ~ 2700 Mq pc -2 
and Ri ~ 190 pc. This is also the approximate radial ex¬ 
tent of the nuclear bar (see Sect. 2). The precise radius at 
which the surface density breaks is Rbreak — 0.75 kpc, and 
the corresponding density is Ebreak — 50 Mg pc^ 2 . This 
means that the nuclear molecular disc extends over four ex¬ 
ponential scale lengths, thus actually more than the main 
molecular disc! 

Our results concerning the nuclear molecular disc of 
NGC 6946 should be compared with those of Regan et al. 
(2001). They showed that the radial profile of CO surface 
brightness decreases steeply for R <1.5 kpc, then increases 
inside a transition layer of radial width A R ss 1 kpc, and fi¬ 
nally decreases again (see their fig. 2). They concluded that 
the CO emission excess in the central kpc looks like a bulge, 
while the decrease at larger radii corresponds to an expo¬ 
nential disc. This is in contrast to our findings. We have 
shown that the CO surface density excess in the central kpc 
is indeed a nuclear exponential disc, and that the transi¬ 
tion from the nuclear to the main disc is sharp: the slope of 
E (R) changes abruptly across the radius Rbreak — 0.75 kpc. 
In addition, the fact that the central ID velocity dispersion 
is very high, ao ~ 50 kms” , does not tell us that the nu¬ 
clear structure is bulge-like, given that the central surface 
density is as large as So « 2500 Mg pc -2 . A real upper limit 
on the central disc scale height can be estimated by neglect¬ 
ing the gravity of the stellar disc, and by assuming that 
the molecular disc is self-gravitating and isothermal along 
the vertical direction: ho <C er 2 /7rGEo ~ 70-80 pc (other as¬ 
sumptions would only modify our estimate by a factor of 
order unity; see, e.g., van der Kruit & Freeman 2011). Since 
the nuclear disc has scale length Ri ~ 190 pc and extends up 
to Rbreak — 0.75 kpC, We filld that ho "C 0.4 Rl W 0.1 Rbreak. 
Thus the central ID velocity dispersion is consistent with a 
nuclear molecular disc of moderate thickness. 

Note that ao ~ 50 kms -1 is comparable to the CO 
values found by Schinnerer et al. (2006) and to the central 
stellar velocity dispersions measured by Engelbracht et al. 
(1996), Ho et al. (2009) and Kormendy et al. (2010). Note 
also that Ro ~ 3 kpc is comparable to the scale length of 
the stellar disc found by Kormendy et al. (2010), and that 
R 1 /R 2 ~ 0.07 is within the range of pseudobulge-to-disc 
scale length ratios observed in spiral galaxies (e.g., Courteau 
et al. 1996; MacArthur et al. 2003; Fisher & Drory 2008). 
As noted by John Kormendy (private communication), these 
results are consistent with the predictions of the general evo¬ 
lution picture reviewed in Kormendy & Kennicutt (2004) 
and Kormendy (2013). 

4.2 Stability of the disc 

Fig. 4 shows the radial profiles of two useful stability diag¬ 
nostics. One is the Toomre-like parameter 




R [kpc] 

Figure 4. Radial profiles of the Q stability parameter (top) and 
characteristic instability wavelength (bottom) for molecular gas 
in NGC 6946. The thick lines are predictions based on the fits to 
k(R), a(R) and E(R) shown in Fig. 3. The thin lines represent 
the approximate radial extent of the nuclear starburst and inner 
bar within bar. 


which measures the stability level of realistically thick gas 
discs (see Romeo & Falstad 2013)0 The other is the char¬ 
acteristic instability wavelength, 

Ac = 2 tt - , (4) 

K 

i.e. the scale at which the disc becomes locally unstable as Q 
drops below unity (see again Romeo & Falstad 2013). Fig. 4 
also shows Q fi t(R) and A c fit(R), which are predictions based 
on the fits to the radial profiles of k, a and E (thick lines). 

The top panel of Fig. 4 illustrates that Q{R) is fairly 
constant for R > 1 kpc, while it varies significantly for 
R < 1 kpc. This is not surprising since the central kpc is 
the region where k(R), <t(R) and E(R) vary most. As we 
move towards the centre, k and a increase by a factor of 
5-7, while E increases by two orders of magnitude (see Fig. 
3). What is surprising is that, in spite of such strong vari¬ 
ations, Q(R) fluctuates around a constant value, deviating 
from it by less than a factor of two! What drives the radial 
variation of Q in the central kpc? We will clarify this point 
in Sect. 5.1. Here instead we focus on the overall flatness 
of Q(R), which is an eloquent example of self-regulation: 

2 We have assumed that the velocity dispersion is isotropic (<r z = 
aj{). This is justified by the fact that interstellar gas can, to first 
approximation, be regarded as a collisional component. 
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the disc of NGC 6946 is driven by processes that keep it at 
a fairly constant stability level, except in the central kpc, 
where the value of Q is modulated by the nuclear starburst 
and inner bar within bar. Current dynamical models of star¬ 
forming galaxies suggest that such self-regulation processes 
are determined by a balance among gravitational instabil¬ 
ity, star formation and accretion (e.g., Bertin & Romeo 1988; 
Romeo 1990; Agertz et al. 2009; Krumholz & Burkert 2010; 
Elmegreen 2011; Cacciato et al. 2012; Forbes et al. 2012, 
2014). According to these models, galactic discs should (i) 
have a constant stability level, and (ii) be close to marginal 
instability. The results illustrated above are clearly consis¬ 
tent with point (i), while the fact that the Q stability pa¬ 
rameter is well above unity (Qat = 4-5) seems to tell us that 
the disc of NGC 6946 is fairly far from marginal instability. 
Is Q really so large? This is a complex issue that we will 
discuss in Sect. 5.2. 

Our results concerning Q{R) seem comparable to those 
of Romeo & Falstad (2013), who showed that the stability 
level of THINGS spirals is, on average, remarkably flat and 
well above unity (see their figs 4 and 5). Note, however, that 
those results cannot be directly compared with ours because 
they concern the THINGS sample as a whole, not each of 
the spirals or in particular NGC 6946. 

The bottom panel of Fig. 4 illustrates that X C (R) in¬ 
creases monotonically for R > 1 kpc, while it has two well- 
defined minima for R < 1 kpc. The location of such minima 
is strongly correlated with the radial extent of the three 
structures in the central kpc: the nuclear starburst (NS), 
the nuclear bar (NB) and the secondary bar (SB). We find 


that: 



-Rminl ^ 

^ 150—200 pc ~ .Rns, Rnb ; 

(5) 

-Rmin2 

- 1 kpc ~ Rsb • 

(6) 


This is surprising because X C (R) is meant to be a diagnos¬ 
tic for characterizing local gravitational instabilities (such as 
those associated with star-forming rings, complexes/clumps 
and starbursts), and not global gravitational instabilities 
(such as those classically associated with the formation of 
spiral and bar structures). Thus Eqs (5) and (6) point out 
Ac (R) as a powerful diagnostic, which accurately predicts 
the sizes of starbursts and bars within bars. In contrast, re¬ 
member that Q{R) is only modulated by such structures. 
What drives the radial variation of A c in the central kpc? 
We will clarify this point in Sect. 5.1. Here instead we go on 
with another important result: 


Ac(-Hminl) 

s .Rns , Rnb ; 

(7) 

Ac(-Hmin2) ^ 

s I Rsb ■ 

(8) 


Note that A C (R) matches the radial extent of both the nu¬ 
clear starburst and the nuclear bar, while it is much shorter 
than the secondary bar. This result indicates that the nu¬ 
clear structures of NGC 6946 are intimately related and were 
most likely generated by the same process of local gravita¬ 
tional instability in the disc, while the strong secondary bar 
has grown significantly after the onset of instability or was 
generated by larger-scale processes. 
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Figure 5. Radial profiles of the Q stability parameter (top) and 
characteristic instability wavelength (bottom) in the inner region 
of NGC 6946. The thick lines show the roles of re, a and E in driv¬ 
ing the radial variations of Q and A c • The thin lines represent the 
approximate radial extent of the nuclear starburst (NS), nuclear 
bar (NB) and secondary bar (SB). 

5 DISCUSSION 

5.1 What drives the radial variations of Q and A c 
in the central kpc? 

In this section, we clarify which physical processes drive the 
radial variations of Q and A c in the central kpc of NGC 
6946. For this purpose, we need to understand in more detail 
how re, a and E contribute to such variations. We can easily 
highlight the contribution of, e.g., re by replacing (a, E) with 
(crfit,Eflt) in Eqs (3) and (4), where ant(R) and Efit(-R) are 
specified in items (ii) and (iii) of Sect. 4.1. Remember, in 
fact, that using these fits for R < 1 kpc suppresses the fun¬ 
damental differences (in a and E) between the central kpc 
and the rest of the disc. The contributions of a and E can 
be highlighted in a similar way. 

Fig. 5 shows the contributions defined above (thick 
lines), together with the original radial profiles of Q and 
A c . The top panel illustrates that re and E play a more sig¬ 
nificant role than a in driving the radial variation of Q: 
re gives rise to the maximum at R ~ 150-200 pc, while E 
contributes to the fluctuations at larger radii. On the other 
hand, a has a significant impact on the magnitude of Q\ it 
helps re to counteract the strongly destabilizing effect of E, 
hence to flatten Q(R)- The bottom panel illustrates that re 
shapes the radial profile of A c more faithfully than that of 
Q. However, even in this case, a has a strong impact on 
the magnitude of A c . In fact, look at the lower curve. Its 
minimum at R ~ 150-200 pc is 4 times deeper than that of 
Ac (R)- This means that ignoring the increase of a towards 
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the centre, i.e. assuming that cr « 10 kms -1 for all R, would 
underpredict the size of the nuclear starburst and bar by as 
much as a factor of 4. 

The bottom line is that radial inflow (as characterized 
by k and E) is the main driver for the radial variations 
of Q and A c in the central kpc. Nevertheless, disc heating 
(as characterized by cr) is of fundamental importance for 
flattening the stability level of the disc and for regulating 
the size of the nuclear structures. 


5.2 Is Q really so large? 


In this section, we discuss the physical factors that can 
push the value of Q down, closer to unity. To estimate the 
magnitude of such effects and impose tighter constraints on 
the stability level of NGC 6946, we need a more general 
Toomre-like parameter. Romeo & Wiegert (2011) introduced 
a simple and accurate approximation for the two-component 
(stars+gas) Q parameter, which takes into account the sta¬ 
bilizing effect of disc thickness and predicts whether the local 
stability level is dominated by stars or gas. Romeo & Fal- 
stad (2013) generalized this approximation to discs made of 
several stellar and/or gaseous components, and to the whole 
range of velocity dispersion anisotropy (a z /ur) observed in 
galactic discs. In the two-component case, the Q stability 
parameter is given by 

W 1 

T 1 Q 1 + T 2 Q 2 

1 W 

T 1 Q 1 + T 2 Q 2 



if T 1 Q 1 > T 2 Q 2 ■ 


if T 2 Q 2 > T\Q\ . 


(9) 


where Qi = K,ai/nGT,i is the Toomre parameter of compo¬ 
nent i, a denotes the radial velocity dispersion, and T) and 
W are given by 


1 + 0.6 (—) 

for 0 < (a z /an)i < 0.5 , 


\urJ i 

(10) 


0.8 + 0.7 (—) 

for 0.5 < ( a z /a R )i < 1 , 


\CR/i 


2cri<72 


(11) 

r i + ^ ' 



This set of equations tells us that the local stability level of 
the disc is dominated by the component with smaller TQ. 
The contribution of the other component is suppressed by 
the W factor. Note, in particular, that if i = 1 is the dom¬ 
inant component and cri differs significantly from 02 , then 
Q « TiQi (as in the one-component case). We are now 
ready to discuss the physical factors that can push the disc 
of NGC 6946 closer to marginal instability. 


(i) Stars. The condition that Q is dominated by stars (*) 
rather than molecular gas (co), T*Q* < T co Q co , can be writ¬ 
ten as 



( 12 ) 


Gerssen & Shapiro Griffin (2012) showed that (o^/ctr)* de¬ 
creases markedly from early- to late-type spirals. For an 
Sc galaxy like NGC 6946, their best-fitting model yieds 
G z /au)* ~ 0.4 (see their fig. 4). Hence T* 1.1. If gas is 
collisional, as is generally assumed, then (a z /an)co ~ 1 and 



0 1 2 3 4 5 

Galactocentric Radius [kpc] 

Figure 6. Radial profiles of the ratios between the surface den¬ 
sities of molecular gas and stars (top) and the ID velocity dis¬ 
persions of molecular and ionized gas (bottom) in NGC 6946. 
The solid lines are the median values of E co /E* (— 0.14) and 
u'co/o'h a (— 0.38) for R < 5 kpc. 


T co ~ 1.5 (see Sect. 4.2). To estimate E*/E co , we retrieve the 
Spitzer IRAC 3.6 fim image observed by the SINGS team 
(Kennicutt et al. 2003), and derive the radial profile of E* 
using the same method as for E co (see Sect. 3.2) and adopt¬ 
ing a fixed if-band mass-to-light ratio, = 0.5 Mq/L q^k 
(Leroy et al. 2008). The top panel of Fig. 6 illustrates that 
E*/E co « 7 for R > 1 kpc, which means that beyond the cen¬ 
tral kpc the stellar and molecular discs have approximately 
the same scale length^ Remembering that er co sw 10 kms" 1 
for R>lkpc (see Fig. 3), the condition that stars domi¬ 
nate over molecular gas finally becomes cr* < 100 kms -1 . 
Since there are no stellar spectra available for NGC 6946, 
we check this condition against the model-based radial pro¬ 
file of <t* proposed by Leroy et al. (2008). We find that 
cr* < 100 kms -1 for i?>2kpc, therefore stars can actually 
dominate the value of Q. But this is not the end of the 
story. To push the value of Q down, from T co Q co « 4-5 to 
T*<3* R3 1, we should have cr* ~ 20 kms -1 , which is true 
only at the edge of the optical disc: R « 11 kpc. Even a 
value of Q « 2 would require cr* « 40-50 kms -1 and oc¬ 
cur at i? « 7 kpc, i.e. still beyond the radial range of our 
analysis. Thus stars alone cannot solve the problem, unless 
the model proposed by Leroy et al. (2008) overestimates the 
real cr*(i?) in NGC 6946. This underlines the importance of 


3 Leroy et al. (2008) found that the stellar disc scale length is 
30% larger than the molecular one, but again their exponential 
fits extend over the central kpc, which is physically distinct from 
the rest of the disc. 
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having observed, rather than model-based, radial profiles of 
a* when analysing the stability of disc galaxies. 

(ii) Atomic gas. Detailed comparative analyses of atomic 
and molecular gas in NGC 6946 show that Shi -C S co for 
R < 5kpc (see fig. 40 of Leroy et al. 2008), while <thi « cr co 
(see fig. 4 of Caldu-Primo et al. 2013). As Thi ~ T co (~ 1.5), 
this implies that ThiQhi T co Q co for R < 5kpc, i.e. atomic 
gas cannot dominate the value of Q within the radial range 
of our analysis. 

(iii) CO-dark molecular gas. There is growing evidence 
that a large fraction of molecular gas in the Milky Way is 
dark, i.e. not traced by CO emission: Mdark « 0.4- 1.6 M co 
(e.g., Paradis et al. 2012; Pineda et al. 2013; Kamenetzky 
et al. 2014; Langer et al. 2014; Smith et al. 2014 and ref¬ 
erences therein; Chen et al. 2015). Can CO-dark molecular 
gas play a significant role in the stability scenario of NGC 
6946? To answer this question, we need at least rough esti¬ 
mates of both Edark and <7d ar k. We estimate Edark from the 
mass fraction found in the Milky Way: Edark ~ 0.4- 1.6 E co - 
The investigations above also suggest that dark molecular 
gas is warmer and more diffuse than the component traced 
by CO emission. So Odark must be bound by <t co (rs cthi) and 
the ID velocity dispersion of ionized gas: a co < Udark < cr Ha . 
We derive the radial profile of ohc, from the Fabry-Perot 
observations of NGC 6946 presented by Fathi et al. (2007), 
following the same method as for o co (see Sect. 3.2). The 
bottom panel of Fig. 6 illustrates that <TH a ~ 2.6 <t co for 
R > 1 kpc, while <thc« approaches <7 C o towards the centre^ 
Hence o co < udark < 2.6 a co . As Td ar k ~ T co {~ 1.5), we find 
that TdarkQdark ~ 0.6-7 T co Q co . This means that CO-dark 
molecular gas can dominate the value of Q and reduce it by 
more than 40%, but only if erdark ~ <t co and Edark ~ 1-6 E co , 
the upper bound observed in the Milky Way (Paradis et al. 
2012 ). 

(iv) Anisotropy of the gas velocity dispersion. As noted 
in Sect. 4.2, interstellar gas is generally regarded as a colli- 
sional component, so that (a z /oR) gas « 1 and T gas ~ 1.5. 
Bournaud et al. (2010) analysed the gas velocity field in 
high-resolution simulations of gas-rich galaxies, and found 
that the velocity dispersion is isotropic only at scales smaller 
than the disc scale height. Agertz et al. (2009) showed that 
at scales of a few 100 pc the gas velocity dispersion has a 
degree of anisotropy similar to the stellar velocity disper¬ 
sion. If this is true for NGC 6946, then ( a z /aR) gas ~ 0.4 
and Tgas ~ 1-1 [see point (i) above], which implies a 30% 
reduction in the value of Q. 

(v) Gas turbulence. Turbulence plays a fundamental role 
in the dynamics and structure of cold interstellar gas (see, 
e.g., Agertz et al. 2009; Hennebelle & Falgarone 2012; 
Falceta-Gongalves et al. 2014; Roy 2015). The most basic 
aspect of gas turbulence is the presence of supersonic mo¬ 
tions. These are usually taken into account by identifying 
(Tgas with the typical ID velocity dispersion of the medium, 
rather than with its thermal sound speed. Another impor¬ 
tant aspect of gas turbulence is the existence of scaling rela¬ 
tions between E ga s, <r g as and the size of the region over which 
such quantities are measured (£). Observations show that 


4 This result agrees qualitatively with that found by Westfall et 
al. (2014), namely that cthc* is on average twice as large as a co in 
nearby spiral galaxies. 


E m - £ 1/3 and (thi ~ i 1 ^ 3 up to scales of 1-10 kpc, whereas 
Eh 2 ~ constant and oh 2 ~ £ 3 ^ 2 up to scales of about 100 pc 
(see, e.g., Hennebelle & Falgarone 2012; Falceta-Gongalves 
et al. 2014; Romeo & Agertz 2014; Roy 2015). Motivated by 
the large observational uncertainties of E gaa (£) and cr gas (£), 
and having in mind near-future applications to high-redshift 
galaxies, Romeo et al. (2010) considered more general scal¬ 
ing relations, E gaa oc £ a and er gas oc £ b , and explored the 
effect of turbulence on the gravitational instability of gas 
discs. They showed that turbulence excites a rich variety of 
stability regimes, several of which have no classical coun¬ 
terpart. See in particular the ‘stability map of turbulence’ 
(fig. 1 of Romeo et al. 2010), which illustrates such stability 
regimes and populates them with observations, simulations 
and models of gas turbulence. Hoffmann & Romeo (2012) 
extended this investigation to two-component discs of stars 
and gas, and analysed the stability of THINGS spirals. They 
showed that gas turbulence alters the condition for star- 
gas decoupling and increases the least stable wavelength, 
but hardly modifies the Q parameter at scales larger than 
about 100 pc. Such predictions have now been confirmed 
and strengthened by the simulation work of Agertz et al. 
(2015). 

(vi) Non-axisymmetric perturbations. Like the original 
Toomre parameter, Q measures the stability of the disc 
against local axisymmetric perturbations, so it assumes that 
kR 1. Here the relevant k is the characteristic instabil¬ 
ity wavenumber, k c = 2n/X c , so the condition kR 1 can 
be written as A c <S 2nR, which is easy to check (see Fig. 
4). While this short-wavelength approximation is satisfied 
by most spiral galaxies (Romeo & Falstad 2013), the as¬ 
sumption of axisymmetric (or tightly wound) perturbations 
is not so general. Local non-axisymmetric stability criteria 
are far more complex than Toomre’s criterion: they depend 
critically on how tightly wound the perturbations are, and 
cannot generally be expressed in terms of a single effective Q 
parameter (e.g., Lau & Bertin 1978; Morozov & Khoperskov 
1986; Bertin et al. 1989; Jog 1992; Lou & Fan 1998; Griv & 
Gedalin 2012). However, there is a general consensus that 
non-axisymmetric perturbations have a destabilizing effect, 
i.e. a disc with 1 < Q < Qcrit is still locally unstable against 
such perturbations. Griv & Gedalin (2012) found that the 
usual estimate Qcrit ~ 2 (see, e.g., Binney & Tremaine 2008) 
is a real upper limit on the critical stability level. This im¬ 
plies that non-axisymmetric perturbations reduce the value 
of Q by a factor l/(3 cr it <S 1/2, i.e. by much less than 50%. 
It is highly non-trivial to constrain the magnitude of this 
effect more tightly than so. The estimates presented in the 
literature (see references above) depend critically not only 
on how open the spiral perturbations are, but also on the 
mathematical treatment of the problem. Tighter constraints 
on the relative importance of non-axisymmetric and axisym¬ 
metric perturbations in spiral galaxies might be found by 
analysing the radial profile of Q for a large galaxy sample, 
and by searching for trends in Q{R) along the Hubble se¬ 
quence. This is however well beyond the scope of the present 
paper. 

The bottom line of this discussion is that stars, CO-dark 
molecular gas, anisotropy of the gas velocity dispersion and 
non-axisymmetric perturbations can all play a significant 
role in the stability scenario of NGC 6946. The disc of NGC 
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6946 can indeed be close to marginal instability. Neverthe¬ 
less, this is most likely the result of all such factors together. 
None of such factors alone can solve the problem. 

Finally, note that the points discussed in this section 
complement, but do not invalidate, the simple and predictive 
stability analysis carried out in Sect. 4.2. In particular, the 
radial profile of Q remains remarkably flat. This is because 
the main effect here is a systematic reduction in the value 
of Q by a factor of order unity. The radial variation of Q 
is less significantly affected, so it is still much weaker than 
that of S, k or a. 


6 CONCLUSIONS 

In this paper, we have provided strong evidence that NGC 
6946 has a double (nuclear+main) molecular disc, with a 
clear transition at R ~ 1 kpc. We have analysed the disc 
structure and stability in detail, using robust statistics as 
well as reliable and predictive diagnostics. Our major con¬ 
clusions are pointed out below. 

• The nuclear molecular disc is exponential, and has cen¬ 
tral surface density Si ~ 2700 Mq pc -2 and scale length 
Ri ~ 190 pc. This is also the approximate radial extent of 
the nuclear bar. The ID velocity dispersion decreases steeply 
from about 50 kms -1 at the galaxy centre to about 10 
kms -1 at the outskirts of the disc, i.e. at a distance of about 
1 kpc. The central ID velocity dispersion is consistent with 
a disc of moderate thickness, given that the central scale 
height is a few tens of parsecs, i.e. an order of magnitude 
smaller than the size of the disc (see Sect. 4.1). 

• The main molecular disc is also exponential, and has 
central surface density S 2 — 66 Mq pc~ 2 and scale length 
R 2 — 2.9 kpc. Note that S 2 is an extrapolated value, since 
the main and the nuclear discs join at a distance of about 1 
kpc. More precisely, the transition in surface density occurs 
at R = .Rbreak — 0.75kpc and S = Ebreak — 50 M e pc~ 2 . 
The ID velocity dispersion is approximately constant up to 
the outermost radius of the analysis presented here (5 kpc), 
and has a value of about 10 kms -1 . 

• These facts imply that the inner region of NGC 6946 is 
physically distinct from the rest of the disc, and that radial 
inflow and disc heating have both played a primary role 
in shaping the central kpc. This is consistent with the fact 
that the dynamics of NGC 6946 is driven by three bars, 
and that the strong secondary bar extends up to R ~ 1 kpc. 
It is also consistent with the classical scenario of secular 
evolution in disc galaxies (see, e.g., Kormendy & Kennicutt 
2004; Kormendy 2013; Sellwood 2014). 

• The radial profile of the Q stability parameter is re¬ 
markably flat. In spite of the fact that E(R), k(R) and 
cr (R) exhibit order-of-magnitude variations, Q(R) fluctuates 
around a constant value, deviating from it by less than a 
factor of two! This is an eloquent example of self-regulation: 
the disc of NGC 6946 is driven by processes that keep it at 
a fairly constant stability level, except in the central kpc, 
where the value of Q is modulated by the nuclear starburst 
and inner bar within bar. The fact that Q is well above 
unity, and on average 4-5, is far more intriguing. We have 
discussed this issue in great detail and shown that the disc 
of NGC 6946 can indeed be close to marginal instability, 
but as a result of four destabilizing factors: stars, CO-dark 


molecular gas, anisotropy of the gas velocity dispersion and 
non-axisymmetric perturbations. 

• The radial profile of the characteristic instability wave¬ 
length, Ac {R), is a powerful diagnostic that accurately pre¬ 
dicts the sizes of starbursts and bars within bars. Such a 
profile indicates that the nuclear starburst and the nuclear 
bar of NGC 6946 are intimately related, and were most likely 
generated by the same process of local gravitational insta¬ 
bility in the disc, while the strong secondary bar has grown 
significantly after the onset of instability or was generated 
by larger-scale processes. 
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